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I Introduction 

The  classical  theory  of  an  underwater  explosion 
bubble  is  based  on  the  assumption  that  the  bubble  re- 
pains  spherical  at  all  times.  However  actual  bubbles 

f • 

do  not  remain  spherical  although  they  may  be  spherical 
initially.  Instead  they  become  flattened  or  even 
kidney-shaped  and  often  break  up.  This  change  of 
shape  is  mainly  due  to  gravity.  The  present  report 
presents  a theory  which  describes  the  change  of  shape 
as  well  as  the  rise  of  the  bubble,  by  taking  account 
of  gravity. 

Several  other  authors,  such  as  H .rring  and  Ward, 
have  also  riven  theories  of  the  change  of  bubble  shape 
^u'-  t''  cravi ty . The  present  wor^  is  intended  to  be 
mere  systematic  and  more  complete  than  any  of  the 
former  theories. 

In  the  present  theory  it  is  assumed  that  the 
water  is  incompressible  and  unbounded,  and  the  pressure 
is  assumed  to  be  uniform  throughout  the  bubble  at  all 
times.  Initially  the  bubble  is  assumed  to  be  spherical. 
Then  if  we  neglect  gravity  we  find  that  the  bubble 
remains  spherical  but  performs  periodic  radial  oscil- 
lations, exactly  as  in  the  classical  theory.  Taking 
account  of  first  o. der  terms  due  tc  gravity,  we  find  that 
the  bubble  remains  spherical  but  also  rises  exactly  as 
predicted  by  Herring’s  formula.  Higher  order  terms 
lead  to  the  change  cf  shape  and  sometimes  to  the  break- 
up of  the  bubble,  and  pIsc  lead  to  a modification  of 
the  Herring  rise  formula.  All  of  these  results  are 
discussed  in  section  VIII  of  this  report. 

» . 

The  assumption  that  the  water  is  incompressible 
can  be  removed,  ad  it  wa3  in  our  previous  report  cn 


p, 

the  spherical  bubble.  However  it  does  not  seem  worth- 
while to  add  the  complication  iOf  compressibility  to  the 
other  difficulties  of  the  present  problem.  Furthermore 
in  section  VIII  wo  have  indicated  how  the  main  effects  of 
compressibility  can  be  obtained  hy  combining  the  present 
results  with  those  of  our  previous  report. 

The  assumption  that  the  water  is  unbounded  can 
alsr  be  removed,  and  the  theory  without  this  assumption 
is  the  subject  of  our  next  report.  In  that  report  it 
will  be  shewn  that  the  results  of  the  present  report 
provide  a first  approximation  to  the  solution  when 
the  bubble  is  not  toe  near  the  boundaries.  Of  course 
further  approximations  will  also  {riven. 

In  section  II  the  problem  is  formulated,  in  section 
III  moving  coordinates  are  introduced  and  in  section  IV 
dimensionless  variables  are  defined.  The  method  of 
solution  is  explained  in  section  V and  the  solution  up 
to  and  including  the  thir^  order  is  considered  in 
section  VI.  A specif. 1 case  is  treated  explicitly  in 
section  VII.  Finally  in  section  VIII  all  the  results  are 
discussed  and  some  general  ccnclusicns  are  drawn. 

II  Formulation 

We  assume  that  an  incompressible  invistid  fluid 
of  infinite  extent  contains  n sras  bubble  within  it. 

The  velocity  u(x,y,z,t)  of  the  fluid  is  assumed  to 
derivable  from  u potential  function  <}>(x,y,z,t)  which 
satisfies  Laplace’s  enuaticn 

(1)  u*  = vi 

(2)  0 

The  pressure  p(x,y,-,t)in  the  wpter  is  then  mi  van  by  the 
Bernoulli  enuation 

(3)  p = m(  t)  -f  f )^j  - fit  z 
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In  (3)  f*  is  the  density  of  water,  g is  the  acceleration 
of  gravity  and  m(t)  is  an  arbitrary  function  of  time. 

We  have  also  assumed  that  the  positive  z axis  points 
vertically  upward. 

The  velocity  ^4  is  assumed  to  vanish  at  infinite 
distance  from  the  bubble 

(!;)  V"4  —*■  0 as  (x2+  y2+  t2)1^2  -*>  co  . 

From  (1)  it  is  clear  that  4 is  defined  up  to  an 
additive  function  of  t.  This  indefiniteness  can  be 
removed  by  specifying  4 or  4t  at  one  point  for  all  t. 
For  this  purpose  we  assume 

(E>)  4 -*  0 as  (x2+  y2+  z2)1/’  -*»co  . 

The  condition  (5)  implies  that  at  infinity,  and 

also  implies  (4),  which  can  therefore  be  omitted. 

The  arbitrary  function  m(t)  in  (3)  can  be  deter- 
mined by  specifying  p at  one  point  for  all  t.  For 
this  purpose  we  assume  that 

2 2 

(6)  p(x,y,o,t)-*  Pq  as  (x  + y ) ->od 

Thus  m(t)  can  be  determined  from  (3,5,6)  which  yield 

(7)  m(t)  = PQ. 

We  call  Pq  the  pressure  at  Infinity  at  z=0.  (It  is 

Just  the  hydrostatic  pressure  at  the  level  z=0, 

namely  Pq=  po+ygZQ,  if  pQ  is  the  atmospheric  pressure 

and  z >o  is  the  depth  of  z=o  below  the  water  surface.] 
o 

The  bubble  surface  is  assumed  to  be  given  by  the 
equation 

(8)  F(x,y,z,t)  = 0 

Since  the  normal  velocity  of  the  bubble  surface  must 
be  the  same  as  thpt  of  the  adjacent  water,  F must 
satisfy  the  kinematic  condition 
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(9)  v'P  • ^ + P = 0 on  F = 0. 

u 

The  pressure  it  within  the  bubble  is  assumed  to 
be  a known  function  of  V,  the  bubble  volume.  ’vc  assume 
that  this  function  is  the  adiabatic  one,  it  = KV~®, 
where  K and  6 are  constants , 6 being  the  adiabatic 
exponent  for  the  gas  within  the  bubble.  Then  because 
pressure  must  be  continuous  across  the  bubble  boundary, 
we  have  the  dynamic  condition  p = KV  ^ or,  using  (3) 

(10)  - f[tt*  ^ (‘Wf!  - pgz  = KV'6  on  F = 0. 

The  mathematical  problem  which  we  consider  is 
that  of  finding  4>  and  F,  satisfying  (2,5,9,10),  given 
F(x,y,z,o)  and  F^(x,y,z,o).  It  is  not  necessary  to 
specify  <|>  initially  since  it  is  determined  by  (2,5,9). 
In  particular,  we  will  assume  that  the  bubble  is  ini- 
tially a sphere  of  raduis  £ and  that  its  initial 

o • 

radipl  velocity  is  a constant  Once  <f>  is  deter- 

mined, u and  p are  given  by  (1,3). 

Ill  Coordinates 

Suppose  that  the  oriein  of  the  x,y,z  coordinate 
system  is  fixed  at  the  center  cf  the  bubble  at  t - 0. 

We  now  introduce  a moving  coordinate  system  by 

the  eouations 

( 11 ) t = x , = y , £ - z -B  ( t ) . 

In  (11)  B(t)  is  a function  to  be  determined  subseauent- 
ly  in  such  a way  that  the  oriein  is  always  at  the 
center  of  gravity  of  the  bubble.  Moving  c^  rdinates 
are  introduced  in  order  that  the  origin  remain  inside 
the  bubble  > s'  long  as  possible. 

Because  of  the  cylindrical  symmetry  of  the 
problem  it  is  advantageous  to  introduce  the  sphorical 
coordinates  r,9,u>.  Obviously  the  solution  will  be 
independent  of  g>,  and  we  will  assume  this  from  the 
outset.  We  define  r,©,w  by 
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(12)  r2=  ^2  + \2+  l ;2  = x2+  rP+(z-B)2  , raO 


cos© 


z-B 


<T7T ?^?v 


0^©<7t 


. I _ X 

tanw  = yr  = - 


0<Cix2tt 


Now  introducing  R(9,t)  we  write  the  eauation  of  the 
bubble  surface  in  the  form 


(13)  F(x,y,z,t ) = r-R(9,t)  = 0 
The  bubble  volume  is  then  given  by 

flirl  V(t)  =—  f11  P'5(e,t)sln8d9 

3 0o 

The  condition  that  the  moving  origin  be  at  the  center 
of  prravity  of  the  bubble  may  bo  written  as 

IT 

(15)  { R4ccs©d©  = 0 

no 

To  express  (2,9,10)  in  terms  of  r and  9,  we 
observe  that  if  is  the  expression  in  terms 

of  for  seme  function  I(x,y,z,t),  i.e.  if 

I(x,y,z,t)  = I(£,>^f£ft),  then  wo  have  the  following 
differentiation  formulae: 


(16)  ix-Tc,  i tT v iz=T^,  it-it4T  «T  4(cob«i  -aineT  ) 

r 


~j'~ ~z  Z.'  _t  t t 

If  we  apply  (16)  to  (9,10)  we  obtain 


(17)  *r— § <J>e  = Rt+fl(cosO+ffi9--RQ)  at  r = R(9,t) 


(18)  K(~  ^ R3sin©d©)"6 


= Pq-  j3  jg(rcos9+B )+<j>£-B( cos* 


sin© 


+©)+  ( V 


- r 


at  r=R(9,t ) . 


6. 


The  problem  is  now  that  of  finding  $,R(©,t) 
and  D(t)  satisfying  (2, 5 ,15 ,17,18)  and  the  initial 
conditions  (19) 


(19)  R(S,o)=A  , R.(©,o)4*  B( o )=B ( o )=0 

U L O 

It  will  prove  convenient  to  employ  the  energy 
equation,  which  is  » direct  consequence  of  the  pre- 
ceding equations,  namely 

IT  71 

(20)  — — C R3alnSde)1'5  + P R3slnSde  - 

6-!  ^ * J0 

Tt 

-TT^C  RP(Rt+B[Cos0+~H®  R0]i(R,©,t)sin©d©  = 


K ,4^3 
6-1  3 ° 


..  P 


^ 


The  terms  on  the  left  side  are  the  internal  energy  of 
the  bubble,  the  potential  energy  and  the  kinetic  energy 
of  the  water,  respectively:  the  terms  on  the  right  are 
the  corresponding  auantities  evaluated  at  t 0,  and 
E is  the  total  energy  which  is  defined  by  this  eauaticn. 


IV  Dimensionless  Variables 

It  is  convenient  to  define  dimensionless  vari- 
ables r,"t  ,T,X,b  by  introducing  a unit  of 

length  L and  a unit  of  time  T.  These  quantities  are 
define^  as  follows: 

~ i»— 6 

( 21 ) r = Lr  B(t)=Lb(t)  k = *o  s 

6-1 

t = T t = L^Po 

4(r  ,©,t ) = lV^F,©,!)  A = La 

T oo 

R(0,t ) = LX(©,t)  A^=  LT”1* 

w O 


6-1  6 
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For  L and  T we  choose  the  definitions 


(22)  L = 


3 ( 6-1  )E 

P7 


1/3 


T = L i 


i P 


T 


\1  o 


Now  equations  ( 2, 5,15,17,18 , 19 ) become  the  following 
set  for  and  b: 


(23)  v24  = 0 , 4 ->  0 


as  r — > 00 


sin© 


(2lf)g(9,t ) = 4r 1 " Xt  " ^cos  6 + 


Xn)  = 0 


at  r = A(©,t ) 


(25)!i(9,t)  ; <Jt  ♦ \ ([1  +(-^)‘;H  i^*h  sineg  + 


0 \2 1 f 1 


V 


+2Xt  -y  f “ 4e  + ^ sin© 3 + X2  - b <y(XP1+  b) 


•tfts! 


\ -6 


A sin©d 9J  -1  = 0 at  r = A(©,t). 


(26)  X(©,o ) = aQ,  At(6,o)  = aQ  , b(o)  = b(o)  = 0 


■K 

(27)  \ 


X'cos©  sin©d©  = 0 


In  (23-27)  the  bars  have  been  omitted  and  two  new 

functions  g(9,t)  and  h(©,t)  have  been  defined  by  (2i|.,25). 

— . 

The  quantity  k is  defined  in  (21).  equation  (25)  is  a 
consequence  of  both  (7)  and  (8),  and  the  symbol  P,  in  it 
denotes  the  Legendre  polynomial  ?^(cos©)  = cos©. 

The  energy  equation  (20)  becomes,  in  dimensionless 
variables  and  omitting  the  bars, 

(28)-  fj-  ^ X2  (Xt+  b(  cos©  + Ae])4(A,e,t)  sin©  d©  + 

"o 

71  It 

+ ^(1  - Cb ) ^ X^sin©  d©  +«-"^6k(^  ^ X^sin©  d©)1“6=oc'^ 


! 
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It  should  be  .noted  that  a and  a are  nc  longer  inde- 

o o 

pendent,  since  they  are  related  by  the  following 
equation,  obtained  by  setting  t = o in  (28). 


(20')|ao  *o 


*1  ♦ ki-3£  =£-3 


Here  &-■  is  the  constant 


(29)  Z = 


6-1  1/2 
6 


V Method  of  Solution 

In  order  to  solve  (23-27)  for  4(r»®»t)»  M©,t), 

and  b(t),  we  make  use  of  the  fact  that  a solution  4 

of  (23)  can  be  expressed  as  a linear  superposition  of 

zonal  harmonics  r“(m+l)p  (cos©)  and  that  \ can  be  ex- 

m 

pressed  as  a linear  superposition  of  Legendre  polynomials 

? (cos©).  The  coefficients  in  these  expansions  are 
m 

functions  of  t.  Next  we  make  the  assumption  that  all 
these  coefficients,  as  well  as  b(t),  can  be  expressed 
as  series (convergent  of  asymptotic)  of  powers  of  the 
parameter  <T.  This  parameter,  defined  in(21),  is  propor- 
tional to  the  acceleration  of  gravity  g,  and  is  equal 
to  the  ratio  of  the  unit  of  length  L to  the  hydrostatic 
head  p0/^g  at  the  initial  bubble  center.  Thus  v/hen 
<T-  0 the  present  solution  is  just  that  of  the  spherical 
bubble  theory  since  (j=  0 is  equivalent  to  g = 0. 

Formally  we  assume 

(30)  4(r,©,t)  = 2 cnm(t)  r"(m+1)P  (cos  9) 

n=o  m=o 

(31)  X(0,t) 

(32)  b(t) 


ifi 

1- 

e* 

QD 

r 

n=o 

m=o 

00 

n 

<5 

XT" 

bn(t> 

n=o 

a (t)  P (cos©) 
nm  m 
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The  problem  is  now  reduced  to  that  cf  determining 

c , a and  b which  arc  all  functions  of  t,  in  such 
run’  nm  n * 

a way  that  (30-32)  satisfy  (24-27),  since  (23)  is  already 
satisfied.  Insertion  of  (31,32)  into  (26)  yields  (dot 
denoting  t derivative) 

(33)  a (o)  = a . a (o)  = a 

oo  ‘ o'  oo'  o 

a (o)  = o,  a (o)  =o  if  n + m > o 

nm  nm 

(34)  bn(°)  = °»  bn(  ° ) ~ °»  n * 0 

We  no’w  insert  (30-32)  into  (24,25)  and  equate  coefficients 
of  -sn  P (cosh),  thus  obtaining  a sequence  of  equations 
involving  the  coefficients  and  their  first  time  deriva- 
tives . 

The  zero  order  terms (n  = c)  correspond  to  the  solu- 
tion obtained  when  C - o,  i.e.  in  the  abo-nco  of  gravity. 
The  equations  for  these  coefficients,  with  the  corres- 
ponding initial  conditions  from  (33)  can  be  solved  at 
once  and  yield 


(35)  c 


om 


a 

om 


o,  for  m > 1 


b 

o 


= o 


The  other  two  zero  order  coefficients  satisfy 

2 • 


(36)  C = - (a  a ) 

^ oo  oo  oo 

2 

(37)  (a2  a )'a~l  - ia 

W OO  00  OO  200 


k 

* 

k 


-36 


oo 


1 


Equation  (37)  is  the  well  known  equation  of  the 
classical  theory  of  spherical  bubbles,  expressed  in  terms 
of  our  units,  end  (36)  is  the  corresponding  expression 
for  the  potential.  The  solution  of  (37)  satisfying  the 
initial  conditions  (33)  is  a periodic  function  of  t 
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which  increases  monotonically  from  a minimum  value  to 
a maximum  value  during  a half  period  end  decreases 
monotonically  to  the  minimum  during  the  next  half  period. 
The  solution  can  bo  written  in  the  form 

(38)  aoo(t)  = +r0  5 

where  the  phase  t is  defined  by  oc-1oc(c'0)  = 


The  function  cx.  ( "f ) is  on  even  periodic  function 
of  C defined  over  the  first  half  period  by  its  inverse 
func  t ion  7T  {cC) 


(39) 


T(  a ) = 


dx 


4 x-3-l-kx-^' 


fit  < a 


The  numbers  X andc*  are  the  minimum  and  maximum 

values  of x t and  arc  the  two  roots  of  the  equation 

(4.0)  l - x3  - kx-3(&-D  = o 


The  parameter  k,  which  determ  nos  a nd'1*-,  is 

defined  in  (21)  and  satisfies  the  inequality 

**  1 

(^1)  0 < k ^ k = — (——=■) 

6-1  'v  0 / 

The  equality  occurs  only  when  aQ  = 0 and  = 

The  second  condition  means  that  the  initial  bubble 

pressure  equals  the  initial  hydrostatic  pressure  at  its 

/\ 

center.  From  (40)  wo  also  have  the  inequality  ( cCis 
defined  in  (29)]. 


(4.2)  0 <£-{k)  < ol<  oe(k ) < 1 


/v 

The  equalities  holds  when  k = k in  which  case  we  have 


/\ 

( 43 ) 0CCC)  =oc=  constant 


I 


11. 

Then  from  (38)  a (t),  which  we  will  henceforth  call 

oo 

a(t),  is  constant  and  equal  to  one.  Thus  the  unit  of 
length  L is  the  equilibrium  bubble  radius,  i.e.  the 
radius  at  which  the  bubble  pressure  would  equal  the 
hydrostatic  pressure  corresponding  to  its  initial 
center,  while  the  total  energy  E would  be  unchanged. 

In  order  to  write  the  equations  for  the  higher 
order  coefficients  (n  * 1)  it  is  corvcnient  to  introduce 
the  functions  g(9,t)  and  fi ( 0 , t ) defined  by 


oo 


(44)  g(e,t)  = g(9,t)  +>__  CT  (b  ?1  [ a .+  aj  anri  + 

n-1  m-o 

+ (m  + 1)  a-(m  + 2)<--  ] P > = £ a-11  £ « P_ 

nn  m °nin  m 


n=o  m=o 


oo 


(45)  fi(0,t)  = h(  9,  t ) -1_  & n(  -36  £ a-36-1  a + 

n=l  no 

+ f a 


A 

k 


7m=o  a2 


+ a a + ] p ) = 

nra  nm  run  m' 


OD  OO 


n=o 


m=o 


h P 
nm  m 


The  functions  g(9,t),  h(6,t)  ccc’urr i/.g  in  (44»45)  arc 

defined  in  (24,25).  The  coefficients  g^,  h^  are  defined 

by  the  above  equations.  It  is  easily  verified  that 

g and  h depend  only  upon  those  coefficients  a„  , 

nm  nm  sm 

ca,  b with  s < n. 
sm*  s 

Now  inserting  (30-32)  into  (24)  we  obtain  for 


n > 0 

e 

nm 


1 am+2 

m+1 


+ 2 — a - z ] 
nm  a nra  &nmJ 


, m 4 1 


(46) 


nl 


1 

7 


a3  [ 


n 


+ a, 


nl 


o a 

2 a anl 


" ®nl^ 


m = 1 


i 


I 


I J 

i 1 


i 


Similarly  inserting  (30-32)  into  (25)  and  making  use  of 


(45,46)  we  finally  obtain  after  some  manipulation 
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(47) 


aS  n + 3aa  + (a  + 36  J a“^6_1)a  = a g.  + 

no  no  k no  °no 


+ 2ag  + h 
°no  no 


(a3bn)-  = - inlJ)  + 2a2  hnl  , 


m = 


m = 


0 

1 


“nm  + 3“nm  + (1*m)  “V.  = °Kn  * (m+2)ienra 


+ (m+l)h  , m>l 

nm 

liquations  (47)  and  the  initial  conditions  (33,34) 

determine  the  a (m4l)  and  b once  a.  and  b,  (k<n)  are 
nm  n km  k 

known.  The  a ..  are  determined  in  terms  of  the  a,  with 
nl  km 

k<n  by  (27).  The  are  given  by  (46). 

Now  all  the  affinal)  arc  zero,  and  bQ=0,  by  (35), 

and  a = a(t)  is  given  by  (38,39).  Thus  we  may  proceed 
oo 

to  find  the  a^  and  b^  and  continue  step  by  step  to 
obtain  all  the  coefficients.  It  would  seem  from  (47) 
that  at  each  step  n>o  we  have  an  infinite  act  of  non- 
homogeneous  ordinary  differential  equations  to  solve. 
However  it  is  shown  in  Appendix  I that  an^,  and  h^ 
are  zero  whenever  m>n  or  m+n  = odd  integer.  Thus  those 
of  equations  (47)  in  which  these  quantities  appear 
become  homogeneous  linear  equations,  and  since  the  initial 
data  are  zero  by  (33,34)  the  solutions  of  these  equations 
are  zero.  Therefore  we  have 

(48)  a ; 0 if  m>n  or  n+m  = odd  integer 

bn  =0  if  n = even  integer 

Thus  from  ( 46 ) 

(49)  = 0 if  ra>n  or  n+m  = odd  integer 

Hence  at  each  step  there  is  only  a finite  set  of  equations 
(exactly  [ ])  to  be  solved. 

As  a consequence  the  expansions  for  4,  and  b. 
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up  to  the  third  power  of  <r",  simplify  tc 

(50)  4 = -a^ar”1  + ^<znr~2?1 ) +<f(C-2Qr~1+C22r~3?2 ) 

+ CT"^C31r  2pi+  33r  ^?3  ^ + * “ 
x = a + cr2{a?0+  a22?2)  + <y3  ( a33p3 ) + •• 

b = <rb^  + + ... 

In  the  expansion  for  X the  terms  and 

were  omitted,  since  by  (27)  f 0,  It  Is  not 

true,  however,  that  a , = 0 for  all  n.  In  fact,  for 
n = 5 «e  h?.vc  a?1  = - -L  *'\2*33  . 

It  i 3 convenient  to  take  advantage  of  the  preceding 

2 

results,  at  least  up  to  terms  in  er  , in  order  to  calculate 
g(0,t)  and  h(©,t ),  which  thon  become 

(51)  g(©,t)  = - 3a  ^2.^a20”  ^ a22^  a ^^^a22P3  + *** 

H(e,t)  = aP1<5'  +(bx  f i b2)^?-2  - ^ b2P2  2 + 

| (5a2Q  + 2a22  + ab1[a“1a22]>1<T3  + 

+ ^ (3a22  " a^1[a’1a22)‘ ) P3^3  + ... 


These  expressions  facilitate  writing  equations  (ij. 7) 
explicitly. 
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VI  Solution  (Jp  to 

For  n = 1 the  only  non-zer^  coefficient  t-c  be 
determined  is  b^.  Since  from  (5l)  wo  have  = o, 

"11  = a»  e<luatien  (47)  yields  for 


(52)  (a^b^ ) = 2c? 

Thus,  us ins  the  initial  data  in  (34) » wo  have 


(53 )  bi  = 2 ( 

"o 

equation  (53)  is  the  Herring  rise  formula  vihick 
accounts  for  the  only  first  order  effect.  The  bubble 
shape  is  unaffected  in  first  order,  but  the  water 
flow  is  modified  since 


d t 


dt 


C = 
11 

For  n = 2 
are  a2Q  and  a22 


dt. 


J 

the^only  coefficients  to  be  determined 
From  (47,51,53)  we  have  the  equations 


(54)  aapp  + 3aap,.  + (a  + 35  a~3b-l)a^  = bi  + ^ bi  # 


a(o)=  ^(o)  0 

rton  uor>  w 


(55)  aCL22  + ^aa22  " aa22 


9 £2  (o)  = ‘(o'/  = 

t D1  ' 22  a 22  G 


In  terns  of  the  solutions  of  these  equations  and 
making  use  of  (I16),  the  only  non- z.oro  coefficients 
<~2in  are  Si ven 

(56)  ^ 20  = **  (a  a20  ^ 

(57)  *~22  = - t a2  (a2aOP) 


22' 


For  n = 3 the  only  coefficients  to  be  determined 
are  and  b^.  From  (47,51,33)  we  have 


I 


(58)  aa^  + 3^^33  - 2aa,~  = 


- 6b^  ( c^2  + ^2?  ^ * 


a^(u ) — c. -*3(0)  — 0 


(59)  (-.3M*  = - 3(a2b1a2C)*+  6r/a20  + | 


b^(c)  = b7('0  --  0 


Then  using  (l^bj 


(6°)  ^ ^j[a3b-  + 3b1a2(a2c  - l a 22)] 


(61)  C. 


- - i(a3 ( i 


bla  :'a223 


The  procodina  equations  o '*  be  .solved  explicitly 
in  the  sooci  1 case  of  the  so-c  lice*  ’"cquilibriusn 
bubble1',  which  is  treat’d  in  the  ne-t  section.  In 
general , however,  it  secr.s  necess  r**  to  integrate 
the  preceding  equations  numerically  in  order  to 
obtain  quantitative  results.  Such  results  ‘.ill  be 
published  soon.  " 

It  is  useful  to  recognize  that  equation  (53) 
for  the  rise  b^(t)  cm  be  gut  in  f •„  fora 

(62)  b1(t)  = S t2  + ? t + 0 

Here  S is  a constant  while*  h and  are  periodic 

functions  of  t havin':  the  saie  period  as  a(t).  This 

form  p omits  us  to  dee  or  nine  the  behavior  of  b^(t) 

for  all  t by  merely  cvalu..ti;v  integrals  for  one  period. 
These  xntc  appear  a . the  oof imtio:i3  of  3,  f 

and  ''  i.hich  arc,  if  t(  is  the  period  cf  a( t ) 

Penny  and  Fr"cc  have  enanineu  f.  solution  of  the  homo- 
feneous  equations  (h7)  for  n = 1 with,  nan-zero  initial 
data  in  order  to  study  the  rrov.’th  of  initial  asvmxotry. 
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(63) 


s = 


I r a3dxi  ' - fto  -*3 


*oJ 


M 

o 


tc  -r  /t 


= 2 \T  { 


i~\ 

_l 


i'3dT  - 


1 ( 0 . 
*o  3. 


-3 


dT 


7 = spt)  - St2  - Pt 


The  periodicity  of  P is  obvious,  while  that  of  can  be 
proved  with  some  manipul  tier:.  hcua^ions  (63)  apply 


only 


r*  • 

= A = a = 0. 


VII  The  Aquilibriun  dubblc 

In  section  V,  foil uWxr-g  c.'.uuti*  n(4l ) it  is  shown 
that  if  k = k then  the  zero  order  bubble  radius  a(t) 
regains  constant  and  equal  to  one.  This  occurs  if 
the  initial  radial  velocity  is  zero  and  if  the  initial 
bubble  pressure  equals  the  pressure  in  the  water  at 
the  depth  of  the  bubble  center.  We  refer  to  this 
case  as  that  of  the  equilibrium  bubble  since  the 
bubble  would  remain  at  rest  in  the  absence  of  gravity. 

oinco  the  radius  u(t)  is  identically  constant 
and  in  fact  equal  to  one,  due  to  our  choice  of  units, 
all  the  equations  of  the  preceding  section  simplify  to 
such  an  extent  that  they  can  be  solved  explicitly. 

3y  solving  them  we  obtain. the  following,  expressions, 
for  X(Q,t),  b(t)  and  ^(r,e,t)  up  to  and  including 
the  third  order  in 


(6 k)  x(e,t)  = i +(crt2)2  jf(r)  - 3A  p2(cosg)j 

3 

(<yt2)  | p^(cose)  + . . 

(65)  b ( t ) =<rt2  - (ert2)  g(r)  + 


• • 
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(66)  4(r,&,t) 


( cozy ) 


}<jV 3 

v5T  / 


+ <?2t^r’^P0(ccs« ) + ( + -"^2  t 
- ^-rr;  sin  ^35  1 1 ' £T ^ r"*^  P.  ( CO30  ) 

(36)V2  / 1 

" r~^  ( ccsO ) 

In  equations  (6*4,65)  7T  = \j3& 't  and  f (7"),  s(7D 
are  given  by 

(67)  f(T)  = - 4r  (1  - cost) 

r-  T4 

(66)  g (tj)  ~ ^77  “ yyn  [2(1  “ ccs(  } -TTsinTj 

These  functions  both  decrease  fairly  slowly  from 
f(o)  = 1/6,  g(o)  = 7/12  to  f(cx>)=  0,  g(co  ) = 9/20  . 

Using  equations  (64)  and  (65)  v;o  have  construct- 
ed graphs  of  the  bubble  profile  at  various  tin.es  from 
t = 0 until  t = 1.C9,  choosing  C = .222.  These 
graphs,  in  Figure  I , 3h.)w  the  rise  of  the  bubble  as 
well  as  its  deformation. 

3incc,  (64, 65)  arc  only  the  first  toms  in  a 
series  solution,  wo  can  expect  them  to  represent  the 
solution  accurately  for  a short  time,  after  which  the 
omitted  terms  become  important.  In  order  to  determine 
when  the  omitted  terms  become  significant,  we  oan 
exa-  inc  the  expressions  in  (64., 65)  to  find  when  they 
fail.  Two  kinds  of  failure  occur.  In  the  first  kind, 
the  origin  of  coordinates,  which  is  located  at  the 
center  of  , gravity  of  the  bubble,  may  cross  the  bubble 
surface  and  entor  the  water.  Then  the  potential 
function  4(r,0,t)  which  has  singularities  at  the  ori- 
gin, will  be  singular  in  the  :ater,  contrary  to  the 


the  solution  was 
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asou  ; tic C'  of  regularity  on  which  the  solution  was 
based.  ihoroforc  the  expressions  (64, &5)  certainly 


beco’ n invalid  as  soon  as  X = 0 for  3 one  values  of 

9 and  t.  If  £->  .222  this  failure  .o  the  first  to 

*■> 

occur,  and  it  occurs  when rr . 2 one.  i'(,r  9 - t. 

Thus  in  this  ease  the  bette  : of  t/.o  tub  bio  bulges 
upward.  until  it  reaches  the  center  of  gravity, 

fhc  sec  «nd  type  of  failure  occurs  because  (65?) 
predicts  that  the  nubble  will  rise  for  awhile  and 
then  fail.  This  is,  of  course,  pir. -si  cully  ucreess.-i- 
ablc,  end  we  therefore  consider  (65)  to  fail  when 


b(t)  = 0,  i.o.  when  the  bubble  stops  rising.  This 
type  of  failure  is  the  first  to  occur  ifeTs  .222, 
and  it  occurs  when  eft  ~.8.  Therefore  vhen  cr  = .222 
both  failures  occur  s i-aultar.eousl",  end  this  occurs 
at  t = l.C-9,  which  exp  la  'ns  why  the  graphs  in  Figure  I 
were  only  constructed  until  this  value  of  t. 

A.s  an  oxa>ivole  of  these  results,  let  us  consider 
a oubble  of  one  inch  radius  initially  located  tea 
inches  below  the  surface.  >i 


U =1  , Z = 10  ) 

o ’ o 

The  pressure  above  the  o'arface  is  at- .ospheric  (p  = 

o 

1 atmosphere)  sc>  that  the  hydrostatic  head  at  the  bubble 

is  33  feet  p lu3  If  inches  or  4O6  inches,  which  is 

40b  tl:-:03  the  initial  radius.  Thus  = .GO246. 

The  unit  of  length  L is  equal  to  the  ir.isxal  radius 

A for  the  equilibrium,  bubble,  30  L = A =1  inch, 
o o 

Then  the  unit  of  tine  f = .0025  seconds . If  we  use 


the  Herring  rise  formula,  which  is  the. first  term  of 

(65),  the  ri..--  b(t)  = C t2,  so  the  bubble  will  reach 

p ZG  . 

the  surface  when  <7t  = — = 1C  or  t = 64  units* 

However  cq  at  ion  (6b)  fa^ls  when  X(ir,t)  = 0 which 
occurs  at  t = 16  units  or  .C4.O  seconds.  The  bubble 
rises  only  .55  inches  before  the  formulae  fail. 
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VI I X D1 3 CU3  3 jtyi ' of  Insults 

•;c  have:  presented  a systc  .atic  -ctkod  for 
determining  the  shape  and  rise  of  the  bubblo  and  the 
potential  function  in  the  vrs.tcr,  iron  ;:hicb  the  pressure 
can  be  computed.  This  method  ic  based  upon  a pov;er 
series  expansion  in  terns  of  a (linen  si  on  leas  para- 
meter or  which  is  the  ratio  of  the  cquilib'  _un  bubble 
radius  to  the  hydrostatic  head  above  the  initial 
bubble  center.  Thus  c~  is  proportional  to  the  accel- 
eration of  gravity  g,  so  our  expansion  car.  also  be 
considered  to  be  in  powers  of  g,  For  any  bubble 
which  docs  not  vent  (i.c.  break  through  the  v.’ate-r- 
air  surface)  during  its  first  expansion,  <S  is  less 
than  one  and  greater  than  zero. 

The  terns  of  zero  order  in  C correspond  to  the 
classical  theory  of  a spherical  bubble.  They  describe 
a solution  in  which  the  bubble  remains  spherical  and 
performs  undated  periodic  radial  oscillations,  while 
its  center  remains  at  rest.  The  equilibrium  bubble 
is  a special  ease  which  does  net  even  oscillate. 

Via  on  terms  of  first  order  are  also  included,  the  only 
modification  is  that  the  center  of  the  bubblo  rises 
according  to  tlio  Herring  rise-  formula, 

\/hen  terms  of  scco.  d order  are  included  the 
shape  of  the  bubble  is  found  to  change  although  the 
rise  formula  is  unmodified.  Tire  changes  in  the  bubble 
arc  of  two  kinds  corresponding  to  the  two  terms 


& a20(t)  and  & aop(t^  ?2(cos«)  in  (50)  respectively. 
The  first  term  corrcsoonds  to  a change  in  volume  of  the 
bubble,  which  is  an  increase  in  the  ease  of  the  equi- 
librium bubble  and  also  in  a numerical  example  which 
we  have  treated,  ouch  an  increase  is  to  be  expected 


in  general  on  physical  grounds,  since  the  bubble  i3 
rising  into  a region  of  lower  pressure.  From  the 


20 


i 


least  initially,  ir.  agreement  with  the  preceding 
considerations*  As  a consequence  •£  the  increase  of 

a^(t)  it  f ollov/3  that  the  tine  between  successive  j 

minima  (or  maxima)  of  the  bubble  volume  is  decreased 

f 

bolov;  that  g iver.  by  the  zero  order  (spherical  bubble) 
tern  alone.  In  addition  the  expansion  phase  as  in- 
creased in  duration,  while  the  contraction  phase  is 
diminished  in  duration. 

The  other  second  order  tern  corresponds  to  a 

i 

flattening  of  the  bubble.  Since  a,^(t)  starts  nc  zero 
and  decreases,  as  can  be  seen  from  the  differential 
cquat ion( 55 ) , this  flattening  is  in  the  vertical 

I 

direction,  i.c.  the  vertical  separation  between  top 
and  bottom  ia  d.ir.i  ished,  while  the  horizontal 

’ t 

separation  between  sides  is  increased.  This  can  bo 
seen  in  the  figures,  and  is  in  afroonent  with  ob- 
servation, a3  are  the  preceding  results. 

When  terns  of  third  order  in  <r  are  included 

I 

the  shaoe  of  the  bubble  is  again  modified  while  the 
rise  formula  is  also  changed.  The  only  tern  of  this 
order  in  the  equation  for  the  bubble  radiu3  is 
C?  *33  ( t)  P^(costt),  which  lcavo3  the  bubble  volume 
unchanged  but  changes  its  shape.  Since  a^(t) 
starts  at  zero  and  increases,  as  one  can  show  from 

rf.'nr.nf  •?  rl  ^ va  f \ tKi  o ’s  r . n r>  f*  olm^r. 

'■'  * w \ ~ - w ^ ^ - — * V 

corresponds  to  pulling  up  the  bottom  raid  pushing  up 

the  too,  so  the  bubble  becomes  kidney-shamed  with  the 

. » - j 

lower  side  c^ncav: . This  behavior  can  me  seen  in  the 

figures  and  i3  al3o  in  agreement  with  observation.  The 
third  order  correction  to  the  rise  formula  can  be  shown 
to  start  at  zero  and  to  become  negative,  from  equation 
(59),  in  agreement  with  the  numerical  oxamplo  and  the 
equilibrium  bubble.  This  correction  i3  also  in  agree- 
ment with  experiment,  since  the  Herring  formula  pre- 
dicts too  large  a rise. 

All  of  the  gravity  effects  described  here  arc 


greater  in  :iGgr.itudc  a: .0  occur  vro  quickly,  the  larger 
th~  value  of  ef,  ’.'m;  <y  is  larger  for  1 argo  explosions, 
so  all  the  ef  fv  cts  described  should.  be  more  pro  -i  icnt 
for  such  explosions . A l.rgo  value  of  <g"  is  also 
obtained  if  the  pressure  above  the  water  surface  is 
reduced.  This  explains  why  these  efforts  have  been 
observed  m reduced  pressure  tanks. 

The  behavior  of  the  bubble  depends  not  only 
upon  c but  else  upon  another  dinens  ionic  a 5 pararictcr 
k.  The  smaller  the  value  of  k the  greater  is  the 
amplitude  ef  oscillation  of  the  bubble,  since  small 
k corresponds  tc  a large  c:qplosion  energy 


The  maximum  value  of  k is  k which  is  attained  for  the 
equilibrium  bubble.  Since  the  gravitational  effects 
depend  mainly  on  C,  result:^  for  the  equilibrium 
bubble  indicate  what  v;ill  tapper,  in  general . There- 
fore we  expect  that  for  any  value  of  k,  equations 
(50)  for  X and  b v;ill  fail  at  a certain  tine,  just  as 
they  do  for  the  equilibrium  bubble,  lc  = k.  Since  in 
genc-ral  the  bubble  becomes  flattened  and  then  kidncy- 
shapod,  wo  expect  the  center  part  of  the  louor  surface 
to  continue  rising,  like  a jet,  until  it  reaches  the 
upper  surface,  thus  converting  the  bubble  into  a 
toru3.  For  the  equilibrium  bubble,  failure  of  the 

p 

equations  occurs  when  Ct  ~.8,  so  wo  may  as  suns  that 
the  real  breakdown  for  any  bubble  will  occur  when 


rl  or  tiS/r 


-1/2 


As  one  oeriod  of  the  bubble  is 


about  3/2  units  for  an  explosion  bubble,  the  breakup 

should  occur  after  about  —---1/2  _ , . , 

30  ' periods.  For  large 

explosions,  with  large  values  ofer,  this  nay  occur 

during  the  first  period. 

If  the  bubble  docs  breakdown  into  a torus,  ■.••hat 

will  be  its  subsequent  behavior?  The  results  of 
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Torso  tika  on  a torus  of  one  fluid  within  another  fluid 
indicate  that  a torus  is  also  unstable-  and  will 
generally  breakup  into  a certain  number  of  pieces 
determined  by  the  dimensions,  etc.  Tlion  each  of  these 
pieces  might  flatten,  become  kidney  shaped  and  break- 
down! into  another  torus,  which  would  again  breakup, 
etc.  This  seems  to  happen  when  a drop  of  one  fluid 
(c.g.  ink)  is  dropped  into  another  (c.g,  water). 

Finally  let  us  consider  hew  compressibility  will 
alter  the  present  results.  Wc  have  attempted  to 
account  for  compressibility  of  the  uatcr,  as  wo  did. 
in  our  previous  report  on  the  spherical  bubble,  by 


using  the  wave  equation  for  the  potential  in  the 
water  rather  than  Laplace’s  equation,^  Uo  then 
expressed  the  solution  of  this  equation  as  an  infinite 
scries  of  multipolos,  but  we  wore  forced  to  use  a 
fixed  origin,  rather  than  a moving  one.  Wc  also 


oiqprcssed  the  bubble  surface  as  a series  in  powers  of  cT . 
and  in  Legendre  polynomials,  and  attempted  to  determine 
the  coefficients.  Wo  found  that  the  terms  of  sere- 
order  in  <5*  were  exactly  the  same  as  the  solution  of 
our  previous  report,  as  we  naturally  expected.  The 
higher  order  coefficients  satisfied  second  order  equa- 
tions similar  to  those  obtained  in  the-  present  report, 
but  the  coefficients  in  those  equations  involved  the 
function  a(t)  of  our  previous  report,  rather  than  the 
incompressible  solution  a(t).  This  was  not  the  only 

uu'w  ^‘vC-.iCCi  bO  L'b  biiu  i;LG3b  l.JjxOI*  bt’.iib  Giitf  • 


he  therefore  expect  that  all  of  the  present  results 


will  still  apply 
the  main  effects 
for  by  replacing 


in  a compressible  fluid,  and  that 
of  compressibility  would  be  accounted 
a(t)  in  all  equations  by  the  compressible 


solution  of  our  previous  report.  In  this  way  the 
main  result  of  that  report-damped  radial  oscillation— 


would  be  incorporated  into  the  present  work. 


.?«. 

"In  an  unpublished  manuscript. 
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Appendix  I 


Theorem:  If  X and  b are  represented  by  the  formal 

series  In  equations  (30-32)  and  If  they  satisfy  the 
equations  (24-27)  then  for  all  n > 0 

(A1 ) g * h = a * c * 0 

Aim  Jim  Jim  urn 

If  m > n or  If  m ♦ n Is  odd,  and 


(A2) 


If  n Is  zero  or  even. 

Furthertaore  and  h^  are  polynomials  In 

’ bk  and  0 < £ < k < n and  are  Independent 

of  all  other  a's  and  b's  except  a... 

oo 

Proof:  The  proof  relies  on  Induction  with  respect  to 

n.  Thus  we  assume  that  the  theorem  is  true  for  all 
n < H and  we  will  show  that  this  Implies  the  theorem 
for  n * N + 1.  It  has  already  been  shown  In  the  text 
that  the  theorem  Is  true  for  n = 0 (see  equations  (35) 
and  (51))  and  therefore  the  theorem  will  follow  for  all  n. 

In  order  to  perform  the  Inductive  step,  we  intro - 

y 

duce  the  set  H£'  of  all  formal  power  series  In  a and  x 
with  coefficients  depending  upon  a parameter  t,  which 
satisfy  the  conditions  u^ft)  * 0 If  n < k,  or  m > n and 
n < £ , or  m ♦ n ’ odd  and  n < £ . The  quantities 
\Xrm( t ) are  the  coefficients  In  such  a formal  series  u 
given  by 


(A3) 


CO  CO 

r *n  z 

n*0  m*0 


-£  £ 

We  also  define  the  subset  C of  power  series  in 

which  the  non-zero  coefficients  u , n < £ , are  polyno- 

nm  “ 

mlals  over  some  set  Co]. 

If  x » cos  0,  then  since  Pn(x)  Is  even  or  odd  ac- 
cording as  n Is  even  or  odd  It  follows  from  the  Induction 
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hypothesis  that 

(A4)  gEHj,  h £ H^. 

If  we  let  [c]  denote  the  set  of  a , c , h and  their 

mo  nm  n 

time  derivatives  with  0 < n < N then  it  also  follows 
from  the  induction  hypothesis  that 

(AS)  J>(r,t)  BhJ,  A £ H*,  bx  £ H^. 

We  first  wish  to  show  that 

(A6)  g £ H^+1,  h £ flj*1. 

In  order  to  do  this  we  first  note  the  following  proper- 
ties of  the  sets  H^: 

0 p \ 

(a)  HjJ'  c §k,  if  JL  » / ’ and  k > k» 

(b)  u £ Hq^  -►  u^a'V0  e H ® n < X 

(c)  u £ 3j^,  v £ Hj^  -►  (Ocu  ♦ pv)  £ , uv  £ fi ^+k 

(d)  u S fi/,  v £ S^\  / + k’  < /’  ♦ k uv  £ 

(e)  u £ H^,  a ^ 0 -*>  (a  ♦ u)n  - nan-1u  £ ft**1, 

( CX  u)n  - a n £ 

(f)  u £ fi/  -►  u.  £ h/ 

O 0 0 

(g)  u £ v £ — ► (1  - x2)uxvx  £ H^+1. 


Now  making  use  of  these  properties  and  the  induction  hypo- 
thesis we  can  prove  (A6). 
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To  prove  (A6)  we  consider  various  parts  of  g and  fi 
separately.  Thus  if  we  consider  the  following  terras  in 
g we  have 


an 


CD 


(A7)  ir(  X,0,t)  ♦ 21  gn  y~  [2  f a^  (m  ♦ l)a 


•(ra+2 ) 


n*l  m*0 
* a2a[ X ~2  ♦ 2a’3( X - a)] 


]P. 


nra  m 


♦ Z Z (»  ♦ 1 )c  (a”(m+2)  - X’(m+2))Pm  £ E^1, 


above  is  in  HN41 
o 


This  is  proved  by  noting  that  since  a ^ 0,  the  first  term 

according  to  (e).  Furthermore  each  term 

in  the  sum  is  in  by  (e),  (d)  and  (a),  and  therefore 

in  H^4^ . Now  by  (c)  the  sum  is  also  in  H^4^.  Thus  the 

statement  (A7)  is  proved.  In  a similar  way  the  remaining 

terms  of  g,  and  all  the  terms  in  H may  be  shown  to  belong 

to  H^41,  but  the  details  will  be  omitted.  Therefore  (A6) 

may  be  considered  proved. 

Now  we  consider  equations  (47)  for  the  a^  and  b 

nm  n 

with  n < N ♦ 1 . Whenever  m>norn+m=  odd,  the  inhomo- 
geneous terms  in  the  equation  for  a_  vanish,  since  these 

nm 

terms  are  proportional  to  h^  and  g^,  which  have  Just 

been  shown  to  vanish.  The  resulting  equations  are  linear 

and  homogeneous , and  by  (33)  the  Initial  data  are  zero* 

Therefore  it  follows  that  the  solutions  a_  =0  for  m > n 

nm 

or  n ♦ m old.  Similarly  from  equation  (47)  for  b, 

find  that  b + a_-  = 0 for  n even, 
n m 


n 


we 


In  order  to  conclude 


that  b 


n 


«•  «• 

4.'/S  V V W*i  n V 


at  prove  that  — 0 for  n 


nl 


even.  Once  this  is  proved.  It  will  also  follow  from  (46) 

that  c__  = 0 for  m > n or  n ♦ m odd  and  the  theorem  will 
nm 

have  been  proved. 

Thus  we  must  show  that  for  N odd,  a^+^  ^ ~ 0.  To 
this  end  we  consider  equation  (27)  and  expand  the  inte- 
grand in  powers  of  cr,  using  (31)  for  X.  The  coefficient 
of  cr**4*  in  this  expansion  is  a sum  of  terms  of  the  form 
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(A8) 


4 


TT 

i*i 


[a 


nimi 


(x)]x  dx. 


In  this  equation  the  variable  x - cos  6 has  been  intro- 
duced. If  all  n^  in  a particular  term  with  a non-vanishing 
integrand  are  < N,  then  since  their  sum  * N + 1 * even, 
the  sum  of  the  corresponding  ra^  is  also  even.  This  fol- 
lows from  the  induction  hypothesis,  since  for  n < N, 

4 


n.  + m.  **  even  if  M 0.  Consequently  TT  (x)  ia 
1 1 nirai  i«l  mi 

an  even  polynomial  in  x and  therefore  the  integral  vanishes 

Hence  the  only  terms  not  vanishing  Identically  are 

those  in  which  one  n,  say  n^  * N + 1 and  the  other  n^  are 


0 since  otherwise  a 


om. 


zero.  Now  if  n^  » 0 then  m^ 

4 x 

Therefore  TT  P (x)  » (x)  and  the  integral  vanishes 

i=l  mi  mJ 

unless  m^  = 1 due  to  the  orthogonality  of  Legendre  poly- 
nomials. Consequently  the  only  non-vanishing  terms  are 
of  the  form 


0. 


(A9) 


(x)x  dx 


2 

3 


aSa 


H+1,1’ 


There  are  four  such  terms  and  thus  the  coefficient  of 
N+l  8 3 

<T  is  ■£  a ajj^  ^ . From  equation  (27)  this  coefficient 
must  vanish  and  therefore 


(A10) 


= °* 


This  completes  the  proof  of  the  theorem. 
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Figure  1 

Rise  and  Change  of  Shape  of  '’Equilibrium”  Bubble 


Profiles  are  shown  for  five  values  of  t from  t * 0 until 
t = 1.89,  at  which  time  breakdown  occurs.  The  center  of 
gravity  is  also  shown  for  each  value  of  t by  a dash  on 
the  z-axis. 
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